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Abstract: We consider QCD at very low temperatures and non-zero quark chemical 
potential from lattice Monte Carlo simulations of the two-color theory in a very small 
spatial volume (the attoscale). In this regime the quark number rises in discrete levels 
in qualitative agreement with what is found analytically at one loop on S 3 x S 1 with 
-R53 <C Aqcd [!]■ The detailed level degeneracy, however, cannot be accounted 
for using weak coupling arguments. At each rise in the quark number there is a 
corresponding spike in the Polyakov line, also in agreement with the perturbative 
results. In addition the quark number susceptibility shows a similar behaviour to 
the Polyakov line and appears to be a good indicator of a confinement- deconfinement 
type of transition. 
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1. Introduction 



The study of QCD at low temperature and non-zero chemical potential /i is compli- 
cated by the presence of the sign problem: setting ji ^ results in a complex action 
which prevents the importance sampling necessary for conventional methods of lat- 
tice simulation. Perturbation theory is also not valid in this regime in the infinite 
volume limit because the coupling strength is large. One method of avoiding the 
sign problem is to study a theory which may closely resemble QCD, specifically the 
two-color theory, for which the action is real even for non-zero chemical potential. 
This is the theory we study in this paper. Another technique which works in the 
low temperature and non-zero chemical potential regime is to compactify the spatial 
volume onto a manifold with size R <C Ag^ D ; this is the limit of the attoworld and 
here it is possible to calculate using perturbation theory. 

In [1] we calculated the phase diagram of QCD in the /J.R-TR plane on S* 3 x 
S 1 , where R is the radius of S* 3 , considering both N c = 3, and also N c = oo for 
which there is a thermodynamic limit and genuine phase transitions can occur. In a 
small volume, quarks occupy well-defined single-particle states of finite degeneracy, 
whose energies are well-separated on the scale defined by T (given by the inverse 
radius of the S* 1 ). In the large N c theory, we found that for RT <C 1 there is a 
rich structure of third-order Gross- Witten-like transitions between "confined" and 
"deconfined" phases (the names refer to the behaviour of the Polyakov line) as the 
chemical potential \i passes one of the quark energy levels. The interpretation is 
that as /i scans past a level, the quarks fill the level and the system becomes de- 
confined in the sense that the Polyakov line gains a non-vanishing expectation value. 
As \i increases further, the system re-enters the confined phase where the Polyakov 
line vanishes. The Polyakov line consequently exhibits a peak with discontinuous 
slope where the phase transitions occur. At finite N c , similar peaks are seen, but 
in this case the curve is smooth since there are no phase transitions away from the 
thermodynamic limit. 

The analysis of [1] was done in the one-loop approximation and it is important 
to consider non-perturbative effects. In this paper we calculate several observables 
using lattice gauge theory simulations of QCD with N c = 2 (also known as Two 
Color QCD or QC2D) formulated on small tori, and at very low temperatures, using 
L 3 x L t = 3 3 x 64 lattices. We will find that the quark number, Polyakov line, 
and quark number susceptibility qualitatively resemble the results from perturbation 
theory on S* 3 x S* 1 [1] but with noticeable quantitative differences, which may result 
from working at larger coupling strength, or may be due to the formulation of the 
theory on a different manifold, an effect which may be possible in the small volume 
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limit in which we work. We supplement further the perturbative results of [1] with 
numerical calculations of the relevant observables on S 3 x S 1 for N c = 2. In neither 
case is there a thermodynamic limit, as both N c and R (L s ) are finite. These findings 
reinforce our observation, made in [1], that deconfining behaviour in gauge theories 
with /i 7^ appears to be associated with a non-zero density of gapless quark states. 

In the next section we introduce the lattice formulation of QCD with gauge 
group SU(2), and in Section ^briefly outline the behaviour of the lattice model on a 
finite system as fi is increased in the non-interacting limit. An important difference 
with the perturbative approach of [1] is that as well as an IR cutoff L s ~ R, there 
is in this case an explicit UV scale associated with the lattice spacing a. Section f| 
presents results from non-perturbative lattice simulations of the interacting theory, 
where possible comparing them with the equivalent quantities calculated on S 3 x S 1 
with N c = 2 and N c = oo. We conclude with a brief discussion. 



2. Formulation and Simulation of the Lattice Model 



We begin by defining the action for Two Color QCD on a hypercubic lattice [2], 
choosing units where the lattice spacing a = 1, and a Euclidean spacetime index 
v = 0, ...,3: 



s = j2 r x M xy \u- ^ - x; tiU > 

x,y,a 







l/Xx ; 



(2-1) 



x,u<\ 



where N c = 2, U u \ is the oriented product of 4 SU(2)-valued link fields U vx around 
the sides of an elementary plaquette in the z/-A plane, i/j, ip are Grassmann-valued 
quark fields located on the lattice sites, whose index a runs over Nf = 2 flavors, and 
fi is the quark chemical potential. The quark matrix utilises the Wilson formulation 
for lattice fermions: 



M xy [U;ii] =S xy - k 



(1 - ju)e^°U ux 5 y)X+ o + (1 + lv )e-^Ut y 8 VjX ^ 



(2.2) 



The parameter = 2N c /g 2 , where g is the bare Yang-Mills coupling, and the hopping 
parameter k is related to the bare quark mass m via 



m 



2k 2k c (P)' 



(2.3) 



In the free field limit — > oo n c = |, but since the action defined via ( |2.2| ) has 
no manifest chiral symmetry, its value is subject to quantum corrections and must 
in general be determined by simulation; chiral symmetry is then only recovered in 
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the limit k — > n c . Chemical potential is introduced via the orthodox prescription of 
treating constant imaginary timelike abelian gauge potential [3]. 



For the system in hand it is possible to calculate quantum corrections non- 
perturbatively using numerical Monte Carlo simulation [2,4]. The simulation pro- 
ceeds via an orthodox hybrid Monte Carlo (HMC) algorithm, which unlike the case 
of QCD is not subject to the notorious Sign Problem for fi ^ 0, since the func- 
tional measure det^M remains real and therefore positive due to the special SU(2) 
property 

KM{^)K~ l = M*(fi) with K = C7572, (2.4) 

where the Pauli matrix T2 acts on color indices. The most important thermodynamic 
observable in the presence of chemical potential is the quark density 



TdlnZ /t 



n q 



The other key observable in this study is the Polyakov line, defined in terms of link 
variables by 

s x C X t=l 1 



3. Free field limit results (f3 = 00) 



In the free field limit U u = 1 the quark density Q2.5|) on a finite lattice can be 
evaluated via a simple mode sum [2]. The results for a 3 3 x 64 system are shown as a 
function of \i for various k in Fig. |l| The numbers are expressed in the form n q /n^ AT , 
where n^ AT = 2N c Nf is the maximum quark density which can be accommodated on 
the lattice as a result of the Pauli principle, essentially one quark of each spin, flavor 
and color per site. The existence of such a "saturated" system is a consequence of 
having an explicit ultra-violet cutoff; there is no corresponding feature in the results 
of [1] which were obtained in the continuum. 

Fig. |l] shows n q increasing monotonically with /1, but with a step-wise behaviour 
rather than the n q oc fi 3 expected in the thermodynamic zero-temperature limit. 
A physical way of understanding this is that for a system with = L t ^> L s , 
the thermal smearing of the Fermi surface associated with degenerate quarks will 
be much smaller than the OIL' 1 ) spacing between plane wave modes, meaning that 
the Fermi surface will be markedly non-spherical [5]. Indeed, close inspection of 
the figure reveals that the total quark occupancy JV = 2N c NfQ on the plateaux 
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Figure 1: Expectation value of quark density as a function of /i for free quarks on a 3 3 x 64 
lattice. 



have respectively Q = 1, 7, 19 and 27, implying the existence of well-separated 
energy levels (on the scale of T) with degeneracies 1, 6, 12 and 8. These are readily 
identified with the plane-wave states available on a 3 3 system, which in units of the 
mode spacing 2vr/L s take the form (0, 0, 0), (±1, 0, 0), (±1, ±1, 0) and (±1, ±1, ±1) 
respectively. The figure also shows that varying k has a small effect on the energies 
of these discrete levels, given by the value of fi at the riser of each step, but none on 
their degeneracy. For the value k = 0.125 corresponding to free massless quarks, the 
onset at which n q rises from zero takes place at fi = 0; this threshold shifts to larger 
fi for k 7^ k c . 

The results for (3 = oo are thus in qualitative accord with the one-loop findings 
of [1], except in that case the equivalent degeneracy of the Lth level is given by 
,jV l = 2N c Nf Yla=i ^ + 1)) with £ = 1, 2, ... as appropriate for angular momentum 
eigenstates on a hypersphere. Another difference is that for massless quarks the 
hypersphere levels q = (£ + |)-R _1 are equally spaced. It is thus reasonable to 
anticipate that the structure of Fig. [I] would remain stable under weakly coupled 
quantum corrections, i.e. with /3 large but finite. 
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4. Results for the Interacting Case 



In this section, we first present numerical results obtained with k = 0.124 at two 
coupling strengths = 24 and (3 = 6, close to the weak-coupling limit. This will 
enable an estimate of how quantum corrections evolve with coupling, but it should be 
noted that both values are considerably weaker than those employed in the studies 
[2,4] which explored significantly larger volumes. The lattice spacing at /3 = 1.9 
was determined via the string tension to be a ~ 0.19fm, implying that L s = 3 
corresponds to a system size considerably smaller than a fermi even at these stronger 
couplings. Using the one-loop beta-function estimate a' /a = exp[— (/?' — f3)/4biN c ] 
with b\ = (llN c —2Nf)/A8iT 2 , we can be confident that the simulations described here 
definitely probe the attoworld. Another important difference is that in previous work 
a diquark source of the form jipip + ji/ji/j (in effect a Majorana mass for the quarks) 
was introduced; this had the effect of mitigating infra-red fluctuations associated 
with the Goldstone modes induced by superfluidity due to diquark condensation, 
thus improving the performance of the HMC algorithm. In the current work we 
judged this would not be needed for such small systems; even so it was found that, 
eg. the number of conjugate gradient iterations required for the HMC acceptance 
step occasionally exceeded 10 3 (see Fig |T0| below) . 




Figure 2: Plot of n q /n q , and quark number susceptibility Xq (rescaled by a factor 
1/25) versus /i for /3 = 24, k = 0.124 on 3 3 x 64. Numbers on the right hand side should be 
multiplied by 2Nj to give the occupancy JV . The meaning of the open symbols is discussed 
in the text 
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Figure 3: Plot of quark number jV , Polyakov line 8? and quark number susceptibility x 
for N c = 2 on S 3 x S 1 . N f = 2, mR = 0, ^ = 50. 



Fig. ^ shows the results of a scan along the /x-axis for /3 = 24. For each value of \i 
2000 HMC trajectories of mean length 0.5 were accumulated. The quark density n q 



is calculated according to eq. (|2.5|) , and we define the quark number susceptibility 
by its derivative with respect to //, 



Xq 



dn q 
dfx 



(4.1) 



In principle in lattice QCD Xq is calculated in terms of correlations of quark bilin- 
ears [6], but in this study since we have a finely-spaced comb of /i- values we adopt the 
more pragmatic route of approximating the derivative in (|4.1|) by a finite difference. 



It is interesting to compare the lattice results with those of a one-loop calculation 
on S 3 x S 1 . Figure |3| shows the same quantities (though normalized differently) 
on S 3 x S 1 and plotted against fiR, where R is the radius of S 3 (for comparison 
< fiL s < 6 in Fig. |2|). The action of QCD on S 3 x S 1 to 1-loop order is given 
as in [1], and was originally derived in [7] for more general matter content. At low 
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temperatures the partition function is given by 
Z(RT) = J [d6] e- s 

/oo 
[d9] exp -]T- [Tr A (P n ) + {-lYN f z f {^mR)e n ^ T Ti F {P n ) 
n=l 

(4.2 

where [dO] = Y[f=i^ii P = diag{e 46>1 , e^ 2 , ...,e tdN }, and Zf(-^^mR) is given by 



z f 



oo 

(^,".fl)=£<4'W" >/T 



(4.3) 



oc 



where m is the quark mass. In Figure ||, the Polyakov line quark number j¥ , and 
quark number susceptibility x-> are an derived from the partition function. These are 
given by 



(4.4) 



Polyakov line: = J [d6] e~ 5 fS^j ' 

fdlnZ\ 

Quark number: jV = T I — - — I , (4.5) 

V 9fi ) 

Quark number susceptibility: y = T I — — ] . (4.6) 

V dfi J 

Note that with this normalisation \ is an extensive quantity. Qualitatively the 
lattice results of Figure || and the perturbative results on S 3 x S 1 in Figure ^| appear 
similar: the fermion number rises in discrete levels, and there is a spike in the 
Polyakov line and the quark number susceptibility at each level transition. However, 
there are noticeable differences in the shape of the curves during the transitions. 
The differences between the lattice and perturbative results could be due to working 
at different interaction strengths. Also, in small volumes effects from considering 
different manifolds should be more apparent. In what follows we discuss each of the 
observables calculated at (3 = 24 on the lattice in detail, comparing with the results 
from S 3 x S 1 where appropriate. 



Before discussing each observable in detail it is interesting to see what happens 
when we increase the interaction strength a little more. Fig. [| shows the same scan of 
the /x-axis at stronger gauge coupling (3 = 6. The level degeneracies corresponding to 
the occupancies on the staircase are identical to those at /3 = 24 indicating that the 
system responds adiabatically in this coupling range. Moreover, the shell energies 
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Figure 4: Plot of n q /n^ AT and versus ji for (3 = 6, k = 0.124 on 3 3 x 64. Results for 
(3 = 24 are also shown as open symbols 



given by the risers are also little changed, showing that dependence of the levels on (3 
is weak. The behaviour of & also suggests a similar relation with n q (n), albeit with 
a much reduced signal which can be attributed to much larger quantum fluctuations 
at this stronger coupling. 

The most striking feature of Figs. ||[| is the qualitative similarity of the relation 
between n q and to that found using weak coupling methods on the hyperspherical 
attoworld of Figure [5] and of the results in [1]. The numerical results reinforce the 
claim that deconfining behaviour is correlated with a partially filled energy level or 
"shell" in the box, implying that deconfinement is associated with a non-zero density 
of gapless states. 



4.1 Quark number density n q 

The weak-coupling results for n q (fj,) in Figure [| exhibit a staircase structure similar 
to that of Fig. [I], but with two important differences. Firstly, the onset value of 
\i ~ 1.2 is considerably larger, and secondly the degeneracies of the discrete levels 
are now 8, 8, 6, 16, 8 and 8, with the understanding that each of these numbers 
should be multiplied by a factor 2Nf (see next paragraph). It should also be noted 
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that over certain /z-ranges n q is double- valued; this arose from simulations where two 
distinct and apparently stable, or at least metastable, states of the system were found 
(in such cases the state with larger n q is shown with open symbols). In some cases 
"tunnelling" events from one state to the other occurred as the simulation proceeded. 
Since both L s and N c are small, and we are accordingly far from any thermodynamic 
limit, there seems to be no obvious criterion for deciding which if either is the "true" 
ground state. 

Before commenting further on the level degeneracies, we wish to highlight an im- 
portant constraint on the interacting theory. The grand canonical partition function 
can be written 

Z(/i)=£z(^exp(^), (4.7) 

where Z is the canonical partition function evaluated in the sector with fixed particle 
number JV = n q V . On one of the steps in Fig. ^ where n g (/i) is approximately 
constant, it must be the case that Z oc exp(^^V), implying that the sum in ( |4.7| ) is 
saturated by one particular value of JV ', and hence that Z oc Z(^¥). Now, using the 
Fourier representation of the 5-function, ([4.7]) may be inverted to yield 



Z{JT) = / d$e- i/r *Z(t$T), (U 

2,71 

i.e. the canonical partition function is the ^Vth Fourier component of the grand 
canonical partition function evaluated with imaginary chemical potential \xi = dT. 
However, it is known [8] that the Z Nc centre symmetry of the pure gauge theory can 
be extended to the case of matter with imaginary chemical potential, resulting in a 
periodicity 

2tt 

Z( l $T) = Z(tT($ + w )). (4.9) 



Eqns. ( |4.7| , f4.9|) together imply Z(.jV) is only non-vanishing for jV jN c integer, or in 
other words that the canonical partition function is only defined for sectors of zero 
iV c -ality. Notice that, while the weak coupling analysis is obviously consistent with 
this constraint, the values of JV observed on the steps of Fig. |2| are, indeed, multiples 
of N c which is a far from trivial check on the lattice analysis. 



In fact, the observed pattern of degeneracies 8, 8, 6, 16, 8, 8 appears impossible 
to explain in terms of single-particle plane wave states, in important contrast to the 
free case Fig. |l|. We conclude that while the staircase behaviour of n q (fi) appears to 
be a universal feature of attosystems, non-perturbative effects continue to play an 
important role in determining the spectrum. 
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Figure 5: Normalized quark number, ^Y/N 2 , on S 3 x 5 1 for the first level transition 



(N c = 2, oo). N = N f = N c , m = 



l 

' TR 



50. 



On S 3 x S 1 the quark number at low temperature and zero quark mass is given 

by [i] 



JV = T 



djj, 



AT r oo N 



=1 i=l 



-ifl i+ =3-(<+I) 



(4.10) 



Here the level structure arises because of the contribution from the Fermi-Dirac 
distribution function 

1 



f(ee) 



(4.11) 



1 _|_ e (e t -iJ.)/T> 

which, at low temperatures, is 1 when fi is larger than an energy level Ei = (£+ 

and zero otherwise. Thus the number of quarks at each level L is a sum over the 

fermion degeneracies 2£(£ + l)N c Nf of all levels with < //: 



oft. 



L 

i=i 



2£(£+l) 



(4.12) 



A plot of the quark number on S 3 x S 1 as a function of /i.R for N c = 2, and N c = oo, 
is given in Figure |5|. In both cases the quark number is symmetric about the mid- 
point of the transition at = £g, in that the level transition starts and finishes at the 
same rate, in contrast to the lattice results of Fig. The N c = oo results also shown 
(from [1]), exhibit quite clearly the two Gross- Witten-type phase transitions as the 
points of non-analyticity where the quark number first turns on and then saturates 
as [iR is increased. 
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4.2 Polyakov line & 



In addition to quark density n q , the other important observable monitored was the 
Polyakov line defined in Eqn. ( |2.6|) . In theories without fundamental matter 3? is an 
order parameter for deconfinement; even in the presence of quarks it can be related 
to the free energy fq of a static fundamental source via & ~ exp(— fq/T), and 
lattice studies of eg. the thermal transition in QCD observe a sharp increase in £P 
around the deconfining temperature. The behaviour of & in systems with \i ^ and 
L t > L s is considerably less studied, although intriguing results have been reported 
in [2,4]. Note that for N c = 2 there is no distinction between and defined 
using the inverse link variables U' in ( |2.(:i| ). 

Fig. |] shows that the Polyakov line exhibits a complicated behaviour: 

in regions where n ? (/i) is constant its value is small, possibly consistent with zero, 
but in regions where n q is changing & > 0, implying deconfinement. In fact it 
exhibits rather sharp maxima over at least six distinct ranges of /i: its numerical 
value &max — 0.4 is in accord with expectations that the theoretical maximum 
value of 1 is renormalised downwards due to quantum fluctuations [9]. When n q 
increases to one of its plateaux over a very short //-interval, the transition from 
2? > to 3? ~ is very sudden; on other occasions when the n q {p) behaviour 
between plateaux is an extended S-shape, &{p) has a more symmetrical peak. This 
seems to be true even in regions where both confining and deconfining states are 
found. Fig. |] focusses on the range 1.23 < /x < 1.33 over which the second shell 
is gradually occupied and L rises from zero to its maximum value. The histogram 
shows that the fluctuations 5^/^ max ~ 0(1) throughout, but that the distribution 
evolves from being symmetrically centred on zero at \i — 1.25 to being highly skewed 
at \x = 1.33. 

The Polyakov lines on S 3 x S 1 as a function of /i.R are shown in Fig. [7| (left) for 
N c = 2 and N c = oo. For N c = 2, = and the Polyakov line is symmetrical 
around each energy level ji = eg. It is important to note that this is not the case for 
N c 2, where ^ and these are individually asymmetrical around each energy 
level, though mirror images of each other. This reflects that fact that in general 
fundamental and antifundamental test charges respond differently to a background 
density of fundamental charge. As remarked above, the lattice results show both 
features which are approximately symmetrical as well as some which are clearly 
skewed, in apparent contradiction with the N c = 2 prediction shown in Fig. 0. 
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Figure 6: Histogram of ^ over a sample of 2000 configurations with (3 = 24 at various /x. 
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Figure 7: Polyakov lines (left) and normalized quark number susceptibility (right) on 
S 3 x S 1 for the first level transition (N c = 2, oo). N = N f = N c , m = 0, ^ = 50. 



4.3 Quark number susceptibility Xq 



To a very good approximation Fig. |] shows that oc % q \ indeed, the quark number 
susceptibility is often taken as an alternative signal of deconfinement in lattice studies 
of the thermal QCD transition [10]. 
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Figure [F] (right) shows the quark number susceptibility from perturbation theory 
on S 3 x S l for N c = 2 and N c — oo. For N c = 2 the quark number susceptibility is 
calculated by numerical integration: 

fdS\ 



X = T- 



—T 



where 



and 



d_ 



h (/ ^ e 



dS_ 



d 2 S 



dS\\ 2 1 /".,„. o fdS\ 2 I 



/V ^ 00 



i=l f=l 

7V C oo 



1 



i(££-/x)-j0i 



8=1 «=1 



1 + eT 



e ^(ei-fJ.)-i8i 



d 2 S 



(4.13) 
(4.14) 

(4.15) 



To obtain the N c = oo result the quark number susceptibility is calculated by 
taking the derivative of the result for the quark number in [1]. It is only non-zero 
while a level transition is taking place; near the L th step the expression is 



2L = T— 

N 2 dp. 



In 



(l+A^)(l + ^j-AA) 



-i 



where Af(p) is determined as in [1] from numerical inversion of 



(4.16) 



(4.17) 



and is given by (|4.12 ). It is clear from eq. (|4.16|) and Figure |7] (right) that for 
N c = oo another derivative of Xg with respect to p will result in discontinuities at 
the beginning and end of each level transition, indicating that these correspond to 
third order transitions as predicted in [1]. 



We note that the level structure of the particle number and corresponding sus- 
ceptibility spikes were also observed in simulations of the non-linear sigma model [11]. 



4.4 Quark— antiquark condensate (ipip) 

Next we discuss another fermionic observable, the quark-ant iquark condensate, which 
requires a non-perturbative calculation. For Wilson lattice fermions (|2.2|) : 

(W) = 4N c Nf - k-^—. (4.18) 
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Figure 8: Plot of versus \i for (3 = 24, k = 0.124 on 3 3 x 64. Open symbols have 

the same meaning as Fig. 0. 



Because ( |2.2| ) has no chiral symmetry, this bilinear has no interpretation as an order 
parameter, but does yield information on the conformal anomaly [2]. The data, 
however, echo the structure of Fig. § (^ip) is approximately constant in /i-ranges 
where there are only filled shells, and changes rapidly in the same regions where n q 
changes rapidly, corresponding to a partially-filled shell. 



On S 3 x S 1 the quark-antiquark condensate is calculated using [1] 



V3 \ dm 
Nfin 



00 N 



T->0 TC 2 B? 



ff 1 „• — 1 



+ 1) 



! Jf,m) 



(4.19) 



which gives zero when m = 0, and a level structure qualitatively similar to Fig. || if 



4.5 Gluonic observables 



It is also of interest to consider gluonic observables. The simplest local gauge- 
invariant gluon observable is the plaquette U v \. In a non-Lorentz invariant system 
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Figure 9: Gluon observables clS Sl function of fi for j3 = 24. (Left): \{U t + U 8 ) (Right): 
□t - □,. 



such as one with fi ^ it is helpful to define 
1 T 



3N r V 



1 T 

3N r V 



(4.20) 

■:').. v, i < — - < — - 

The normalisations ensure \3 s ,t — > 1 in the /3 — > oo limit. We then consider in 
Fig. |] both the average plaquette + and the difference — D s , which is 
proportional to the gluon energy density 

e g = 3Zp{n t - □,), (4.21) 

where Z is a renormalisation factor calculable in perturbation theory [12]. 



The average plaquette has a value extremely close to unity, as befits such a 
weak coupling. It has a non-monotonic variation with //, qualitatively similar to 
the behaviour found on larger systems in [2,4]. The decrease at large /i is readily 
understood as a consequence of saturation n g /n^ AT —> 1: in this regime screening 
due to virtual quark - antiquark pairs is suppressed due to Pauli blocking, and hence 
many gluonic observables revert to their values in the quenched theory, which has in 
effect a larger lattice spacing and hence a larger departure from the free-field value. 
At /i w 1.7, where the average plaquette dips below its value at /i = 0, the ratio 
n q/ n q AT ~ 0.7, to be contrasted with the corresponding value ~ 0.018 observed at 
f3 = 1.9 [4]. 

The difference Dj — D s also shows non-monotonic behaviour; its negative value 
makes it hard to interpret as a physical energy density, and probably arises as an 
artifact of the lattice aspect ratio L t ^> L s . Both plots hint at a finer structure, such 
as a dip at fi ~ 1.65. Curiously this does not appear to match any interesting region 
of Fig. 0, except in the sense that both confining and non-confining solutions appear 
to be stable here. 
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Figure 10: Number of congrad iterations required in the acceptance step of the HMC 
algorithm versus jj, for j3 = 24, k = 0.124 on 3 3 x 64. Blue and green points correspond to 
the open symbols of Fig. ^. 



4.6 Conjugate gradient iterations 



Finally, Fig. [10] plots the number of iterations of the conjugate gradient algorithm 
required to invert M to some pre-specified accuracy (in this case such that the average 
norm of the residual vector should not exceed 10~ 9 per site, spin and color) as a 
function of \i. This quantity carries physical information since it is roughly inversely 
proportional to |A m j n | 2 , where \ m %n is the complex eigenvalue of M closest to the 
origin. Fig. [10] suggests that the smallest eigenvalue evolves smoothly so long as 
n q (fj) is smooth, but jumps sharply at the same points where n q (fi) is discontinuous. 
Note that the eigenvalue is always smallest in deconfined regions where \qi & > 0> 
again consistent with the picture of a non-zero density of gapless modes. Interestingly, 
this pattern is also respected by the alternative solutions corresponding to the open 
symbols of Fig. || suggesting that at these points the simulation sometimes evolves 
smoothly with fi, and at other times jumps to a new ground state with significantly 
smaller X m i n ; further studies on systems with varying L t (correpsonding to varying 
T) might clarify the situation here. 
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5. Summary 



This paper has followed up the study of QCD at non-zero chemical potential on 
spaces defined by a scale R <C ^qcd initiated in Ref. [1] by presenting Monte Carlo 
results for the case N c = 2 which is simulable using orthodox lattice gauge theory 
techniques. We have shown that the two principal qualitative features of the pertur- 
bative calculation on the hypersphere, namely that the quark density n q is a step- wise 
function of //, and that the Polyakov line is significantly different from zero only in 
the //-ranges where n q is rising, persist on the three-torus {S 1 ) 3 x S 1 even once non- 
perturbative effects are correctly included. The staircase form of n q (fj) is indicative 
of a series of widely separated sets of physical states akin to the shell structure of, say, 
a nucleus. Moreover we have identified a further relation, oc \ q = dn q /dfi, which 
seems to be respected equally well in both approaches. These results lend support 
to the interpretation of \ q as an alternative indicator of confinement/deconfinement 
in systems with fundamental matter, and hint at a relation between deconfinement 
and a non-vanishing density of gapless states (in other words, the existence of a 
"conduction band"). 

However, the two approaches differ in their detailed predictions for the shell 
degeneracies. Comparison of simulation results at two different (3 values suggests 
the system evolves adiabatically with coupling, but it has not proved possible to 
interpret the levels in terms of the single-quark states underpinning the perturbative 
approach, and a full explanation of the structure revealed in Fig. |2] is still missing. 

In addition, identification of the correct ground state of the system as /x is varied 
was troublesome due to the apparent existence of more than one "solution" stable 
under HMC evolution (in fact, this proved to be a much more serious problem in 
pilot studies on 2 3 x 64, and persuaded us to switch to a system where L s was not a 
multiple of N c ). In conventional simulation campaigns such ambiguities are normally 
resolved by taking the thermodynamic limit, but here L s and N c are both finite. It 
would, of course, be interesting to try a simultaneous tuning of L s and (3, keeping 
L s a fixed, but this requires resources considerably beyond what we have been able 
to expend. In any case, more general questions about how both thermodynamic and 
zero temperature limits are approached would also be interesting to explore. 
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